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ABSTRACT 

We present the results of a comprehensive re-analysis of the images of a virtually 
complete sample of 28 powerful SCR radio galaxies with redshifts 0.6 < z < 1.8 from 
the HST archive. Using a two-dimensional modelling technique we have derived scale- 
lengths and absolute magnitudes for a total of 16 3CR galaxies with a median redshift 
of z = 0.8. Our results confirm the basic conclusions of Best, Longair & Rottgering 
(1997, 1998) in that we also find z = 1 3CR galaxies to be massive, well-evolved ellip- 
ticals, whose infrared emission is dominated by starlight. However, we in fact find that 
the scalelength distribution of SCR galaxies at z ~ 1 is completely indistinguishable 
from that derived for their low-redshift counterparts from our own recently-completed 
HST study of AGN hosts at z ~ 0.2. There is thus no evidence that SCR radio galaxies 
at z ~ 1 are dynamically different from SCR galaxies at low redshift. Moreover, for a 
10-object sub-sample we have determined the galaxy parameters with sufficient accu- 
racy to demonstrate, for the first time, that the z ~ 1 SCR galaxies follow a Kormendy 
relation which is indistinguishable from that displayed by low-redshift ellipticals if one 
allows for purely passive evolution. The implied rather modest level of passive evo- 
lution since z ~ 1 is consistent with that predicted from spectrophotometric models 
provided one assumes a high formation redshift (z > 4) within a low-density Universe. 
We conclude that there is no convincing evidence for significant dynamical evolution 
among SCR galaxies in the redshift interval < z < 1, and that simple passive evo- 
lution remains an acceptable interpretation of the K — z relation for powerful radio 
galaxies. 
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1 INTRODUCTION 

It has been known since the early 1980's that the hosts of 
powerful radio galaxies display a tight relation between the 
T^-band magnitude and redshift (Lilly & Longair 1984). In 
recent years it has been shown that, at least out to 2 = 1, 
essentially the same relation is followed by both the less pow- 
erful Parkes Selected Regions (PSR) and 6C radio galaxies 
(Dunlop et al. 1989; Eales et al. 1997) and by X-ray selected 
brightest cluster galaxies (Collins & Mann 1997) . When the 
A'-band magnitudes of the SCR galaxies are corrected for the 
expected effects of passive evolution of their stellar popula- 
tions they appear to represent rather good standard candles 
from the present day out to redshift one and greater. Con- 
sequently, for many years it was widely accepted that the 
tightness of the K — z relation for the SCR radio galaxies 
could be most naturally explained by the most powerful ra- 
dio galaxies having a rather well-defined mass, being formed 
at z 2> 1, and evolving basically passively thereafter. 

However, this long-held view has recently been chal- 
lenged in a series of papers by Best, Longair & Rottgering 



(1997,1998), hereafter BLR. Between 1994-1996 BLR un- 
dertook an extensive study of a virtually complete sam- 
ple of 28 powerful FRII SCR galaxies in the redshift range 
0.6 < z < 1.8 (see Fig |). BLR obtained V and I-band HST 
images, J and K-band UKIRT images, and complimentary 
radio observations at 8.4 GHz with the VLA. 

BLR made use of the four broad-band images they ob- 
tained for each of their sources to perform spectral synthesis 
fitting. The four broad-band fluxes where fitted by a simple 
two-component model consisting of an old stellar population 
and a power-law contribution which represented any possible 
aligned component. The stellar population SED's were con- 
structed from the models of Bruzual & Chariot (199S), and 
assumed a 1-Gyr burst of star formation at Zfor = 10, with 
the stars evolving passively from that point onwards to the 
redshift of the particular SCR galaxy. The assumed form of 
any possible nuclear/aligned emission was a power-law of the 
form fi, OL which was used as a compromise between 

the likely contributions of quasar/starburst/dust scattered 
emission. The results of this fitting procedure (BLR 1998) 
show that the broad-band fluxes from these SCR galaxies 
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Figure 1. The redshift distribution of the eomplete 79 objeet 
FRII SCR subsample (Laing et al. 1983) . The shaded objects are 
those included in the Best et al. sample. Reproduced from Best 
et al. 1997. 



are fully consistent with that expected from an old passively 
evolving stellar population, as required in the traditional in- 
terpretation of the K — z relation. 

However, the results of the determination of the ellipti- 
cal galaxy haJf-light radii (t g) performed by BLR, combined 
with a summary of the literature concerning the environ- 
ments of SCR galaxies at high and low redshift, led BLR to 
challenge the traditional view of the K—z relation. Using the 
method of fitting dc Vaucouleurs templates to azimuthally 
averaged luminosity profiles, BLR produced best-fit scale- 
lengths for 19 of their 28 objects (excluding 3C41 and 3C22 
which have a large nuclear contribution (Leyshon & Eales 
1998) not dealt with in the BLR fitting scheme). The re- 
sults of this profile fitting revealed a mean scalelength of 
re = 14.7 ± 1.3 kpc (^o = l,i?o = 50). As pointed out by 
BLR this is significantly larger than the re = 8.2 ± 1 kpc 
found to be typical of low redshift ellipticals by Schombert 
(1987), and only around a factor of two smaller than the av- 
erage scalelength of \aw-z brightest cluster galaxies (BCG), 
Ve = 32.7 ± 1.1 kpc, found by the same author. Combining 
this result with the La r^'^ relation found by Kormendy 
(1977), and confirmed by our HST AGN host galaxy study 
(Dunlop et al 2000), BLR argue that even at 2 = 1 the hosts 
of 3CR radio galaxies are highly evolved massive systems, 
possibly larger by a factor of two than their \ow-z counter- 
parts. 

Combining the apparent similarity of the z — 1 3CR 
galaxies to BCG at the same redshift, with results from the 
literature concerning the environments of high redshift radio 
galaxies (see BLR 1998 for a review), BLR concluded that 
the passive appearance of the K — z relation is a cosmic con- 
spiracy. The general picture given by existing studies is that 
the environments of 2 ~ 1 3CR's are consistent with mod- 
erately rich clusters; galaxy-galaxy cross-correlation func- 
tion measures are consistent with Abell class 0, or richer 
(Hill & Lilly 1991), while multi-colour imaging shows many 
companion objects consistent with passively evolving coeval 
galaxies (Dickinson 1997). As pointed out by BLR, the same 
is not found in existing studies of the environments of low 



redshift 3CR's (Prestage and Peacock 1988), where it has 
been reported that these objects reside in fairly low-density 
environments and have smaller characteristic sizes and lu- 
minosities than their high-redshift counterparts (Lilly and 
Prestage 1987). 

On the basis of the apparent redshift dependence of 
3CR cluster environments BLR argued that the 3CR galax- 
ies at 2 = and 2 = 1 cannot be linked by a simple "closed- 
box" passive evolution scheme. In the picture favoured by 
BLR, the K — z relation can be explained by the two 3CR 
galaxy samples having different dynamical evolutionary his- 
tories, but being observed at the points in their respective 
histories where they both contain a characteristic mass of 
stars; a few times 10^^ M©. The 2 ~ 1 3CR galaxies would 
reach this mass first due to their location in a larger density 
peak producing strong merging between z = 3 ^ 1. The 
3CR galaxies at \aw-z would reach the characteristic mass 
later through merging within their weaker cluster environ- 
ment between 2 = 1 and the present day, as is expected 
in many hierarchical galaxy formation models (Kauffmann 
1999). The SCR galaxies at 2 = 1 would obviously also grow 
through merger activity in this redshift range and would be- 
come the radio dormant BCG we see today at low redshift, 
radio activity having ceased presumably due to the lack of 
available gas to feed their central black-hole. Although this 
picture is both consistent and appealing, it is dependent 
on the acceptance of environmental studies which now date 
back 15 — 20 years, together with the finding that SCR galax- 
ies at 2 ~ 1 have sizes comparable to BCG at that epoch, 
something which is not found to be true for \aw-z SCR galax- 
ies. 

It was a careful examination of the methods used by 
BLR in determining their scalelengths which provided the 
original motivation for this paper. As mentioned above, 
scalelengths were presented by BLR for 19 of their 28 ob- 
ject sample (the images of the excluded 9 objects having 
poor signal-to-noise or sufi'ering from extreme aligned emis- 
sion). However, only 6 of these 19 scalelengths where ac- 
tually determined from the high resolution HST imaging. 
The reason for this is that the one-dimensional profile fitting 
technique employed was unable to cope with the masking of 
significant contamination by companion objects or aligned 
emission. The majority of the scalelength information pre- 
sented by BLR was in fact derived from the analysis of their 
Jf-band imaging data obtained at UKIRT. As reported by 
BLR, this imaging was obtained in 1" seeing conditions with 
typically 54 minutes of on-source integration which, consid- 
ering the redshift range of the objects, leaves a relatively 
small amount of data which has both good signal-to-noise 
and is free from substantial seeing effects. The use of an 
analytical gausian function to represent what has proven to 
be the extremely compUcated IRCAMS PSF (McLure et al 
1999b) casts further doubt on the quoted scalelength values. 
Experience gained from the comparison of our modelling re- 
sults of 2 ~ 0.2 AGN host galaxies (McLure et al 1999a, 
McLure et al 1999b) with those previously obtained from 
typically 1" ic'-band imaging (Dunlop et al. 1993, Taylor et 
al. 1996) has demonstrated to us that the twin benefits of 
high resolution and a temporally stable PSF tend to combine 
to revise the best-fitting galaxy scalelength downwards. The 
main reason for this systematic trend appears to be that the 
high spatial resolution provided by HST allows the reliable 
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Figure 2. Histogram siiowing tiie redsiiift distribution of the 
virtually complete Best et al. sample. Objects shaded grey and 
black are objects which have been succesfuUy modelled during 
this re-analysis. Objects shaded in black are members of the lO-o 
bject sub-sample. 

identification of separate companion objects and their sub- 
sequent masking from the modeUing process, preventing the 
inclusion of extra flux from biasing the scalelength determi- 
nation to higher values. This effect is obviously only going 
to be strengthened at redshifts of z ^ 1 where the angular 
separation between companion objects and the target can 
easily be of the same order as the typical seeing experienced 
during the BLR UKIRT observations. 

Given these considerations it was felt worthwhile to 
re-examine the publically available HST data to determine 
whether the use of the full two-dimensional analysis utilised 
in the 2 ~ 0.2 quasar host galaxy programme would pro- 
duce significantly different scalelength figures from those 
published by BLR. 



2 DATA REDUCTION 

All the HST images taken of the SCR galaxies in the 28- 
object BLR sample were obtained from the HST archive 
facUityQ A detailed list of redshifts, filters and exposure 
times can be found in Table 1 of BLR (1997). An investi- 
gation of the shorter-wavelength exposures of each object 
(mainly F555W, F622W) confirmed that they were either 
of insufficient signal-to-noise to be useful, or dominated by 
emission aligned with the radio axis. Of the 28 objects in the 
sample, preliminary analysis of the /-band images revealed 
a total of 16 which could be successfully modelled with our 
two-dimensional technique. The observational parameters of 
these 16 objects are detailed in Table |l] and illustrated as 
part of the full sample in Fig |^ . 

For each object there are two CR-SPLIT 7-band expo- 
sures of unequal length available, together giving a typical 
exposure time of ~ 1800 seconds. The initial processing of 
the images, flat-fielding and bias removal, was carried out 
by the standard HST pipeline. The two exposures were then 
combined using the iraf task CRREJ, which successfully re- 
moves cosmic ray events using a sophisticated sigma-clipping 
algorithm. The next step in the reduction process was the 

* http: / / archive.stsci.edu 



fitting of a plane to the image with the SCR galaxy masked 
out, to accurately determine the sky background while al- 
lowing for any residual flat-fielding gradients that may have 
been present. The final step was the production of a two- 
dimensional mask for each source which eliminated any com- 
panion objects, or aligned emission, from the model-fitting 
process, as well as any regions of the image that could have 
been biased by scattered light from nearby bright stars. 

In the vast majority of cases the high resolution pro- 
vided by the HST allowed the identification, and subsequent 
masking, of companion objects and diffuse aligned emission 
in a straightforward fashion. For a small number of sources 
(SC252 & SC239) the question of which areas of the images 
should be masked was more ambiguous, due to the HST im- 
ages revealing double nuclei which are not apparent in near- 
infrared images (BLR 1997). As mentioned by BLR, the 
detection of two nuclear components in these sources sug- 
gests either that they are the products of recent mergers, or 
alternatively, contain a central dust lane. Given the undis- 
turbed appearance of the K-band images of these sources 
(BLR 1997), it was decided that no attempt would be made 
to mask the central regions of the HST images before mod- 
elling. 

None of the objects comprising the 10-object sub- 
sample (see Section |^ appear to be highly disturbed in the 
HST images, and consequently return model fits which are 
largely insensitive to the details of aligned-emission masking. 
The advantage of being able to mask areas of the HST im- 
ages, while still running the model fitting in two-dimensions, 
has allowed us to take greater advantage of the HST data 
than was possible with the one-dimensional analysis tech- 
niques employed by BLR. 



2.1 Empirical PSF Determination 

Simple radial surface-brightness plots for the 16 objects cho- 
sen for modelling reveals that there is sufficient signal-to- 
noise to allow fitting to a typical radius of ~ 4" . It is there- 
fore necessary to have a point spread function (PSF) to con- 
volve our model galaxies with which has good signal-to-noise 
out to at least this radius. Due to the synthetic PSF's pro- 
duced by the tinytim software package (Krist 1995) being 
unable to reproduce the WFPC2's scattered light halo out- 
side a radius of ~ 1.5", it is clear that an empirical PSF is 
required for the modelling of these data. A glance at Table 
|l] reveals that all of the exposures utilised here were imaged 
through either the F785LP or F814W filters. This required 
the acquisition of two relatively deep PSF's imaged with the 
correct filter/chip combination, which were also located close 
to the average chip position of the SCR galaxies, in order 
to avoid the noticeable positional variation of the WFPC2 
PSF. An interrogation of the HST PSF search tool |T| pro- 
duced disappointing results, with all available PSF's either 
too faint or too small in angular extent. 

Fortunately, two suitable stars were present on the ex- 
posures of SC41 (F785LP) and SC239 (F814W). Due to the 
need for sufficient depth in the PSF wings, both of these 
stars had saturated cores in even the shortest exposures. 

t http:/ /www. stsci.edu/instruments/wfpc2/ 
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Source 


z 


Filter 


Exp Time/s 


III 

'^1/2/ 


»^l/2/kpc 


Ml/2 




J IT 

^nucf ^host 


3C22 


0.938 


F814W 


1400 


0.7 


5.9 


22, 


,3 


19.8 


0.306 


3C34* 


0.690 


F785LP 


1700 


3.0 


23.8 


24, 


.2 


18.8 


0.002 


3C41* 


0.795 


F785LP 


1700 


0.8 


6.3 


22, 


.2 


19.7 


0.170 


3C49* 


0.621 


F814W 


1400 


2.0 


15.0 


23, 


.0 


18.9 


0.000 


3C65 


1.176 


F814W 


1760 


1.2 


9.9 


24, 


.5 


20.9 


0.088 


3C217 


0.897 


F814W 


1700 


0.7 


5.8 


22, 


.0 


20.1 


0.007 


3C226* 


0.820 


F785LP 


1700 


1.0 


8.5 


22, 


.6 


19.4 


0.050 


3c239 


1.781 


F814W 


2200 


1.3 


10.8 


23, 


.5 


19.7 


0.041 


3C247* 


0.749 


F814W 


2400 


3.1 


25.2 


24, 


.1 


18.6 


0.025 


3c252 


1.105 


F814W 


1700 


0.9 


7.3 


23, 


.2 


20.5 


0.000 


3C277.2* 


0.766 


F814W 


2400 


0.9 


7.6 


22, 


.5 


19.7 


0.005 


3C289 


0.967 


F814W 


1800 


2.1 


18.0 


24, 


.9 


20.0 


0.063 


3C337* 


0.635 


F814W 


1400 


0.8 


6.4 


22, 


,5 


19.8 


0.055 


3C340' 


0.775 


F785LP 


1700 


0.5 


4.3 


21, 


,3 


19.7 


0.000 


3C352' 


0.806 


F814W 


1800 


2.2 


18.1 


24, 


,7 


19.6 


0.024 


3C441* 


0.708 


F785LP 


1700 


1.4 


10.9 


22, 


.7 


19.1 


0.023 



Table 1. The results of the two-dimensional modelling of the 16-objects from the Best et al. z ~ I 3CR sample which did not suffer 
from excessive aligned emission. Listed in column 5 are the fitted scalelengths in arcseconds, with the corresponding figure in kpc for 
Oo = li ^0 = 50 listed in column 6. The surface-brightness at the fitted scalelength in Cousins I-magnitudes arcsec"^ is given in column 
7. Column 8 lists the integrated apparent I— band magnitude of the best-fit host galaxy with column 9 giving the ratio of the integrated 
luminosity of the best-fit nuclear component and host galaxy. Objects labelled with a * are members of the 10-object sub-sample. 



Source 


z 


A»l/2 




3C34 


0.690 


21.7 


-24.9 


3C41 


0.795 


19.4 


-24.5 


3C49 


0.621 


20.5 


-24.5 


3C226 


0.820 


19.7 


-24.9 


3C247 


0.749 


21.2 


-25.4 


3C277.2 


0.766 


19.5 


-24.4 


3C337 


0.635 


19.9 


-23.7 


3C340 


0.775 


18.5 


-24.4 


3C352 


0.806 


21.1 


-24.7 


3C441 


0.708 


20.1 


-24.7 



Table 2. The absolute magnitudes and characteristic surface- 
brightness of the 10-object sub-sample {Hq = 50, Qo = 1) 



To overcome this problem a modified version of the PSF re- 
sampling technique described in McLure et al 1999b was im- 
plemented. Both PSF's were only saturated within a radius 
of < 0.3" of their core, well inside the radius where tinytim 
can accurately reproduce the empirical PSF. Making use of 
this, the two PSF's had their core replaced with the equiv- 
alent TINYTIM model, the relative scaling being determined 
by matching the fiux in an annulus between 0.4" < r < 0.7". 
Another advantage of this approach is that due to tinytim 's 
ability to produce model PSF's at up to 50 times oversam- 
pling, sub-pixel centring of the PSF can be matched to that 
of the galaxies to an accuracy of < 0.005". Although not as 
crucial as with the quasar host work it is still important to 
overcome the severe undersampling of WFPC2, especially 
as several of the SCR galaxies modelled proved to have a 
substantial point-source contribution. 



3 MODELLING 

The modelling of these objects proceeded in an identical 
fashion to that of the z ~ 0.2 AGN sample studied by 
McLure et al. (1999a) and Dunlop et al. (1999). Each object 
was modelled with both a standard Freeman disc template 
and de Vaucouleurs r^''* profile, with no a priori assump- 
tions being made about position angle, axial ratio or any 
possible central point-source contribution. Prior to the mod- 
elling of each source a maximum radius, outside of which pix- 
els would be excluded for insufficient signal-to-noise, was de- 
termined by an examination of the one-dimensional surface- 
brightness profile. The median value for this outer radius 
was 4", corresponding to a physical radius of 33 kpc at the 
median redshift of z — 0.8 {Ho = 50, Qo = !)• This is reason- 
ably well matched to the outer radius of 12" (50 kpc) used 
for each of the radio galaxies in the z — 0.2 AGN sample, en- 
suring that the same physical region of the galaxy light dis- 
tribution is under investigation in both redshift regimes. The 
modelling technique strongly preferred a r^'''* host galaxy for 
all 16 objects analysed. The quality of the fits was good, with 
no significant residual flux. 



4 RESULTS 

In this section the scalelength, absolute luminosity and Ko- 
rmendy relation results from our two-dimensional modelling 
of the z ~ 0.8 objects are presented. Due to the significant ef- 
fect that choice of cosmology can have on angular diameter, 
cosmological dimming and look-back time over this redshift 
range, we give these results for a range of possible cosmolo- 
gies. Four different representative scenarios are considered 
featuring two values of Ho (50,70), together with both open 
(fio ~ 0.1) and flat (fio = 1) geometry. 

The results from the modelling of all 16 HST objects are 
presented in Table |^. The main conclusions reached from 
this modelling work are based largely on the results for a 



© 0000 RAS, MNRAS 000, 000-000 



The Evolution of SCR Radio Galaxies from z=l 5 




o 

a 

tn 
>^ 
C/2 



7000 8000 9000 
A/Angstroms 

Figure 3. The F814W filter profile including the system response 
and quantum efliciency of WFPC2. Also shown is the location of 
the 4000A break of an elliptical galaxy spectrum observed at a 
redshift of 2 = 0.9 



CO 



O 
CO 





Ho 


median / kpc 


< re > / kpc 


1 


50 


9.7 


12.6 ±2.3 


1 


70 


8.2 


9.3 ±1.6 


0.1 


50 


11.5 


14.8 ± 2.7 


0.1 


70 


8.2 


10.6 ± 1.9 



Table 3. The scalelength results from the two-dimensional mod- 
elling of the z ~ 0.8 sub-sample. Columns one and two detail the 
choice of cosmology. Column three gives the median scalelength 
values for the 10-objects in the sample. Column four gives the 
corresponding mean values together with the standard error. 
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10-object sub-sample. The sub-sample consists of the low- 
redshift end of the BLR sample and is illustrated in Fig ^ 
The remaining 6 objects (z > 0.9) have been excluded from 
the following analysis due to the incursion of the 4000A 
break into the HST F814W filter. The throughput of the 
F814W filter is shown in Fig |^ complete with the system 
response and CCD quantum efficiency. As can be seen from 
this Figure, for objects with redshifts z > 0.9 this filter will 
bridge the 4000A break, meaning that it can no longer be 
assumed that the detected flux has originated from the dom- 
inant old stellar population. Given that the alignment effect 
in radio galaxies is stronger in the rest-frame UV than at 
longer wavelengths it is to be expected that these objects 
are contaminated by a signiflcant aligned component. The 
10-object sub-sample has the advantage of having a tight 
redshift distribution with a mean of 0.74 ± 0.07 and a me- 
dian of 0.76, which is well matched by the distribution of the 
10 radio galaxies studied at low-z, with its mean redshift of 
0.20 ± 0.04 and median of 0.20. This allows a direct com- 
parison between the two samples without the added compli- 
cation of allowing for significant evolutionary effects within 
the samples themselves. 

4.1 Scalelengths 

The results of our determination of galaxy half-light radius 
(re) for the 10-objects in the z ~ 0.8 sub-sample are pre- 



Figure 4. A comparison of the distribution of scalelengths de- 
rived from our modelling of the z ~ 0.8 3CR galaxies with that ob- 
tained by Best et al. Also shown is the distribution of scalelengths 
obtained for the 10 radio galaxies and 10 radio-loud quasars 
from our HST z ~ 0.2 AGN host-galaxy programme (Dunlop 
et al. 1999; McLure et al. 1999a, b). All three histograms assume 
Ho = 50, = 1 



sented in Table |3|. Two features of this table are immediately 
obvious. Firstly, it is clear that regardless of cosmology the 
derived half-light radii are remarkably consistent, with both 
the mean and median values agreeing to within 4 kpc. Sec- 
ondly, the mean scalelength is systematically larger than the 
median in all four cosmologies, consistent with the sugges- 
tion from Fig ^ that the distribution of scalelengths from the 
full 16-object sample, and indeed our 20 radio-loud AGN at 
z ~ 0.2, has a substantial tail toward high values. 

Given that the scalelength results obtained by BLR are 
quoted assumming Ho = 50, Qq = 1.0 it is straightforward to 
investigate what differences exist between their scalelength 
determinations and those presented here. Due to the fact 
that most of the BLR scalelength information is derived 
from their K-band observations rather than the HST im- 
ages, it is not the case that they have published a scalelength 
value for each of the objects re-modelled in this paper. In 
order to perform a comparison between the two sets of re- 
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Ho 


z 


Median Mi 


< Mi > 


1.0 


50 


0.2 


-24.52 


-23.45 ±0.16 


1.0 


50 


0.8 


-24.60 


-24.62 ±0.13 


1.0 


70 


0.2 


-23.69 


-23.72 ±0.16 


1.0 


70 


0.8 


-23.87 


-23.89 ± 0.13 


0.1 


50 


0.2 


-24.54 


-24.55 ± 0.16 


0.1 


50 


0.8 


-24.98 


-25.07 ±0.18 


0.1 


70 


0.2 


-23.80 


-23.81 ± 0.16 


0.1 


70 


0.8 


-24.25 


-24.24 ±0.13 



Table 4. The absolute /-band magnitudes of the two 10-object 
radio galaxy sub-samples. Column four gives the median figures 
with column five listing the corresponding mean figures complete 
with standard error. 

suits it has been necessary to simply base the BLR figures 
on the 8-objects from our 10-object sub-sample for which 
they have derived a scalelength value. 

Assuming Ho = 50, f2o = 1-0 the mean scalelength of 
the 10-object sample is = 12.61 ±2.26 kpc with a median 
of 9.70 kpc. The corresponding BLR-derived values for 8 of 
these 10 objects are = 15.15 ± 2.71 kpc with a median of 
14.95 kpc. Given the small number statistics that are avail- 
able, and the inherent difficulty in constraining galaxy scale- 
lengths, the median is the more robust measure of the typical 
scalelength. Using the median it can be seen from these fig- 
ures that the one- dimensional analysis technique employed 
by BLR has systematically overestimated the characteristic 
size of the SCR galaxies by « 50%. 

In order to check that the use of these two cut-down 
samples was not overtly biasing the results, a comparison 
was also performed between the full 16-object sample and 
the 12-objects from this sample for which there is also a 
BLR-derived scalelength figure available. The mean scale- 
length of the 16-object sample is — 11.48 ± 1.59 kpc with 
a median of 9.20 kpc. The corresponding BLR-derived val- 
ues for 12 of these objects are = 15.66 ± 1.91 kpc with a 
median of 15.45 kpc. It can be seen from this that the use of 
the expanded samples strengthens the conclusion that the 
BLR scalelengths are overestimated, with the median BLR 
scalelength being ~ 70% greater than the two-dimensional 
modelling results presented here. 

4.2 Absolute Magnitudes 

The absolute Cousins /-band magnitudes for the z ~ 0.8 
and z ~ 0.2 sub-samples are listed in Table ^ for the four 
different cosmologies. The values shown are calculated from 
integrating the best-fit de Vaucouleurs profile to infinite ra- 
dius in order to be consistent with the published results for 
the low-2; AGN hosts (McLure et al 1999a, Dunlop et al 
1999). As was pointed out in Section |^ the objects compris- 
ing the sub-sample were imaged through two separate HST 
filters; F814W and F785LR The F814W filter is a mem- 
ber of the standard HST filter set and closely mimics the 
Cousins /-band filter. As a result, it was decided to convert 
all of the integrated apparent magnitudes to their equivalent 
Cousins /-band value. The 5 objects from the sub-sample 
imaged with the F814W filter had no correction applied to 
them since the similarity of the F814W and Cousins /-band 



filters is such that for elliptical galaxies the difference in 
magnitude is expected to be less than 0.05 in all cases. This 
level of photometric accuracy is substantially greater than is 
possible in the face of the uncertainties in determining the 
host magnitudes. 

The throughput for the F785LP filter is significantly 
different from that of the Cousins /-band filter. The conver- 
sion from magnitudes obtained through this filter to Cousins 
/-band magnitudes is therefore more complicated. The orig- 
inal strategy to overcome this problem was to make use of 
the one object from the 16-object sample (3C239) for which 
there are exposures in both filters available with comparable 
signal-to-noise. However, this object has the highest redshift, 
z — 1.781, of all of the 28 objects in the BLR sample and 
subsequently the two images were too faint to get a reliable 
conversion. In addition to this, 3C239 has a highly distorted 
morphology in the HST images (BLR 1997), making it un- 
clear whether a conversion factor obtained from this object 
would be applicable to the mostly undistorted objects in the 
sub-sample. In light of this, the method used to convert the 
F785LP fluxes for each host to the equivalent F814W fig- 
ure, and hence Cousins I-band magnitude, was to predict 
the count rates Robject (e~^s~'^pixel~^) of a source of ap- 
parent visual magnitude V in the two filters using the filter 
ratios of Holtzman et al 1995. 

The cosmology-independent k-corrections for each ob- 
ject, produced by the blueward shifting of the /-band fil- 
ter along the galaxy spectrum with increasing redshift, were 
calculated from the figures presented for a burst-elliptical 
galaxy by Rocca-Volmerange & Guiderdoni (1988). Given 
that the predictions produced by different spectral synthe- 
sis codes can differ significantly (Chariot et al 1996) it was 
considered worthwhile to make an independent check of the 
validity of these k-corrections. Due to the restriction of this 
analysis to the 10-object sub-sample, for which the images 
sample galaxy light longward of the 4000A break, it is pos- 
sible to estimate the necessary cosmological k-corrections 
by modelling the galaxy spectrum as a power-law of the 
form f^ai'~°', where a is the spectral index. The value of 
a appropriate for radio galaxies imaged in the /-band was 
estimated to be a ~ 2 from a typical old (> 12Gyr) ellip- 
tical galaxy spectrum. Alternatively, the mean I — J colour 
of the sub-sample objects can be used to estimate the spec- 
tral index, a process which also results in a figure of a ~ 2. 
Both of these estimates are in good agreement with the value 
of a = 1.82 needed to reproduce the k-corrections of Rocca- 
Volmerange and Guiderdoni (1988). The question of whether 
the absolute magnitudes presented in Table |^ are substan- 
tially different from the results previously obtained for the 
z = 0.2 radio galaxies is addressed within Section |^. 

4.3 The Kormendy Relation 

Given that the hosts of the 10 radio galaxies within the sub- 
sample are well fitted by a standard de Vaucouleurs galaxy 
template, it is interesting to see whether the parameters ob- 
tained from these fits produce a Kormendy relation compa- 
rable to that followed by both low-z inactive ellipticals (and 
the hosts of the z = 0.2 AGN, see Section ^). The ^1/2 val- 
ues required to construct the Kormendy relation have been 
corrected for the cosmological dimming of surface-brightness 
according to: 
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/^l/2 


= 3.51±o.36 log ''1/2 + 16.53±0.38 


no 


= 1.0 


^^0 


= 50 


/^l/2 


= 3.51±o.36 log r-i/2 + 17.05±o.33 


no 


= 1.0 


Ho 


= 70 


/^l/2 


= 3.51±o.38 logr-1/2 + 16.28±o.42 


no 


= 0.1 


Ho 


= 50 


Ml/2 


= 3.51±o.38 log r-1/2 + 16.80±o.37 


no 


= 0.1 


Ho 


= 70 



Table 5. The best-fitting Kormendy relations for the 10 z ~ 0.8 
3CR galaxies under four different choices of cosmology. 



no 


Ho 


median r'e / kpc 


< re > / kpc 


1 


50 


11.3 


13.8 ±2.2 


1 


70 


8.1 


9.8 ±1.6 


0.1 


50 


11.9 


14.4 ±2.3 


0.1 


70 


8.5 


10.3 ± 1.6 



Table 6. The scalelength results from the two-dimensional mod- 
elling of the 2 ~ 0.2 sub-sample. Columns one and two detail the 
choice of cosmology. Column three gives the median scalelength 
values for the 10-objects in the sample. Column four gives the 
corresponding mean values together with the standard error. 



h 
h 



(1) 



(l+2l)3+" 

where a value of a = 1.8 has been assumed. The resulting 
~ ''1/2 relations for the four cosmologies are given in 
Table |5| It is clear from this that the choice of cosmology 
makes no significant difference to the slope of the Kormendy 
relation (as expected given the small redshift range of the 
objects), and that in all cases this slope is consistent with 
that of ~ 3 displayed by inactive low-z ellipticals (Kormendy 
1977). 

It is worth noting that the ri/2 and /ii/2 parameters 
from the one- dimensional modelling of BLR failed to pro- 
duce a /ii/2 — '"1/2 relation which was consistent with the 
expected slope of ~ 3. The modelling results of BLR had 
the 2 ~ 1 SCR galaxies lying along a constant luminosity 
slope of 5, exactly as is expected when the galaxy luminosi- 
ties have been well determined but the scalelengths have not 
been constrained (Abraham et al. 1992). The fact that the 
relation presented above is consistent with the Kormendy 
relation, within the errors, can therefore be taken as fur- 
ther evidence that our two-dimensional modelling has been 
much more successful in constraining the scalelengths of the 
sub-sample objects. The possibility of using the Kormendy 
relations derived above to test for the eff ects of passive or 
dynamical evolution is explored in section 5.3. 



5 COMPARISON WITH LOW-REDSHIFT 
RADIO GALAXIES 

5.1 Scalelengths 

The derived median and mean scalelengths for the 10-source 
z ~ 0.2 radio galaxy sample from our low- 2; AGN host galaxy 
study (McLure et al. 1999a, Dunlop et al. 1999) are pre- 
sented in Table H. A comparison of these results with those 



presented for the z ~ 0.8 sub-sample in Table |^ shows the 
two groups of galaxies to have very similar characteristic 
scalelengths. The poorest agreement between the two sets 
of results is for the Q = 1.0, Ho = 50 cosmology. How- 
ever, even in this case the figures are consistent to within 
the errors. In the other three cosmologies the median and 
mean scalelengths for the two samples are basically identical, 
differing by < 0.5 kpc in all cases. The similarity between 
the scalelength distributions of the two sub-samples is con- 
firmed by an application of the Kolmogorov-Smirnov (KS) 
test. For the two Q,o = 1-0 cosmologies the KS test returns 
a probability of 0.68 that the two samples are drawn from 
the same underlying distribution. This conclusion is even 
stronger in the Qo ~ 0.1 cosmologies where the KS test re- 
turns a probability of 0.97 that the two distributions are 
the same. Therefore, we find no evidence that the galaxies 
comprising the z ~ 0.8 sub-sample are systematically larger 
than their z ~ 0.2 counterparts. 



5.2 Absolute Luminosity 

In order to allow a direct comparison of the characteristic 
luminosities of the low- and high-redshift radio galaxy sub- 
samples, and thereby look for any evidence of significant 
merger activity, it is first necessary to correct for the differ- 
ent filters and the expected passive evolution of the stellar 
populations. To achieve this it was decided to convert the 
z ~ 0.2 i?-band magnitudes to their /-band equivalents, 
and then to predict the brightening of the stellar popula- 
tion between z — 0.2 and z = 0.8 due to passive evolution 
alone. The present-day colour of an elliptical galaxy formed 
at high redshift {z > 3) has been taken as R — I — 0.7 
(Fukugita et al. 1995). The corrections to be made for the 
effects of passive evolution have been calculated using the 
synthetic galaxy spectral models of Jimenez et al. (1996) 
using our chosen set of four possible cosmologies. Several 
different galaxy formation redshifts were considered to inves- 
tigate their effects on the resulting correction. A comparison 
between the predicted amount of /-band passive evolution 
between z — 0.2 and z = 0.8 and the difference in median 
absolute luminosity of the two samples is presented in Table 



I 



It is immediately clear from these results that the dif- 
ference in absolute /-band luminosity between the z ~ 0.2 
and z ~ 0.8 samples is inconsistent with the amount of 
passive evolution predicted by the stellar synthesis models 
within the two Einstein-de Sitter cosmologies. This is still 
true even in the fio = 1 model which allows present-day el- 
lipticals to be as old as possible (z/or ~ 10, //o = 50); this 
model still requires 2.5 times more magnitudes of passive 
evolution than is seen in the data. 

Adoption of either of the open cosmologies results in a 
measured difference in absolute luminosity of ~ 0.5 mags 
between the two galaxy samples (columns 3 & 4). As can 
be seen from Table ^, this luminosity difference is formally 
consistent with the expected passive evolution in both open 
cosmologies, for all alternative formation redshifts. However, 
it is interesting to note from Table |^ that the Ho — 70, flo ~ 
0.1 cosmology clearly provides the best match between the 
data and the spectrophotometric model predictions. In this 
cosmological picture the model predictions are perfectly con- 
sistent with the data for all star formation redshifts of 2 > 4, 
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0.18±0.16 


0.52 


0.60 


0.70 


1 


70 


0.17 


0.18±0.16 


0.63 


0.60 


0.56 


0.1 


50 


0.52 


0.44±0.16 


0.26 


0.24 


0.46 


0.1 


70 


0.43 


0.45±0.16 


0.38 


0.51 


0.64 



Table 7. The results of the tests to determine the influence of different cosmology and formation redshift upon the amount of passive 
evolution expected between z = 0.8 and z = 0.2. Columns one and two list the cosmological parameters for each test. Column three 
lists the amount of evolution required to reconcile the median of the two absolute magnitude distributions. Column four gives the 
corresponding mean figures complete with standard errors. Columns 5-7 list the amount of passive evolution between z = 0.2 and z = 0.8 
predicted by our spectrophotometric modelling for three galaxy formation redshifts. (The apparently anomalous trend displayed by the 
f2o = Ij Ho = 70, 2/or = 3 model is due to the fact that even at z = 0.2 a galaxy is only 6 Gyr old, and remains relatively bright at I.) 





= 2.86±o.34 log ri/2 + 17.66±o.38 


Ho 


= 1.0 


^^0 


= 50 


/^l/2 


= 2.86±o.34 logri/2 + 18.08±o.33 
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= 70 
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= 50 
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Table 8. The best-fitting Kormendy relations for the 10 2; ~ 0.2 
radio galaxies under four different choices of cosmology. 
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Ho 


AM{-la) 


AM 


AAf(+lcr) 


1.0 


50 


-0.24 


0.43 


0.56 


1.0 


70 


-0.24 


0.44 


0.58 


0.1 


50 


0.18 


0.60 


0.84 


0.1 


70 


0.19 


0.60 


0.84 



Table 9. The amount of luminosity evolution predicted by forcing 
both galaxy sub-samples to lie on a Kormendy relation with slope 
3.2. Column four lists the implied vertical shift required to over- 
lay the Kormendy relations of the two samples using the best-fit 
scalelengths of Table ^ converted to the appropriate cosmology. 
Columns 3 and 5 list the vertical shift implied by shifting the best- 
fit scalelengths left(— Icr) and right{-flcr) until the scalelength 
distributions of the two sub-samples differ at the la level. 



in good agreement with recent discoveries of old stellar pop- 
ulation ellipticals at redshifts of z = 1.5 — > 2 (Dunlop et al. 
1996, Spinrad et al. 1997, Stiavelli et al 1999 ). 



5.3 Kormendy Relation 

If the two samples of radio galaxies can truly be linked by a 
single population of passively evolving ellipticals it is to be 
expected that they should follow Kormendy relations which 
are identical except for a simple vertical shift in surface- 
brightness. After ma king the appropriate surface-brightness 
corrections (Section O) and R —> J-band filter transfor- 
mation (Section 5.2) the least-squares fit to the Kormendy 
relations formed by the z ~ 0.2 radio galaxies are shown in 
Table |. 



It can be seen from this that, as for the 2 ~ 0.8 galaxies, 
the choice of cosmology makes little difference to the slope 
of the Kormendy relation, although the normalization does 
change significantly. It is also apparent that the low-z ob- 
jects appear to follow a substantially fiatter relation, with 
a slope of 2.9 instead of 3.5, although clearly both values 
are consistent with the expected slope ~ 3 (due to the fairly 
substantial formal error in the fitted slope, which is pre- 
dominantly a result of the small sample size and lack of dy- 
namic range in radius). To facilitate a fair comparison of the 
high- and low-z radio galaxy Kormendy relations, the least- 
squares fitting was repeated with an enforced intermediate 
slope of 3.20 . If it is indeed the case that these two galaxy 
populations can be linked by passive evolution alone then 
the vertical shift required to reconcile the Kormendy rela- 
tions should be in good agreement with both that required 
to match the absolute magnitudes, and the passive evolu- 
tion predictions of the stellar synthesis modelling. To ex- 
plore what range of vertical magnitude shifts are allowed by 
the data, we performed the least-squares fitting (with fixed 
slope=3.20) not only with the best estimate of re for each 
galaxy, but with the scalelengths of the z ~ 0.8 sub-sample 
shortened and lengthened such that the KS test showed its 
scalelength distribution to differ from that of the z ~ 0.2 
sub-sample at the la level. The results of this process are 
presented in Table |^. A comparison of the figures from Table 
P with those of Table ^ shows that in the D,o = 1.0 models 
considered here it is only possible to reconcile the luminos- 
ity evolution predicted by the Kormendy relations with that 
predicted by the spectral modelling and the absolute mag- 
nitude distributions, separately. In these cosmologies it is 
impossible to force both alternative measures of the amount 
of luminosity evolution, and the predictions of the spectral 
modelling, into agreement. In contrast, the ilcr predictions 
from the Kormendy relations in the Slo ~ 0.1 cosmologies 
comfortably bracket the offsets given in Table |^, thus leav- 
ing both alternative measures of the luminosity evolution in 
excellent agreement. 

In theory, the Kormendy relations for the low and high-z 
sub-samples presented here offer an opportunity to constrain 
both cosmology and the prevalence of merger activity. How- 
ever, in practise the effects of these are very closely coupled. 
For example, our results can be reproduced either by pure 
passive evolution since z ~ 0.8 in an open Universe, or by 
modest growth (~ 20% growth in scalelength & luminos- 
ity) since z ~ 0.8 in an Einstein-de Sitter Universe. Despite 
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Logio(^i/2 / kP'=) Log^o(r^/2 / kpc) 

Figure 5. Shown in the top-left panel is the Kormendy relation followed by the z ~ 0.8 sub-sample which has a best-fit slope Of 3.5. 
In the top-right panel the z ~ 0.2 radio galaxies (crosses) have been added along with their best-fit relation of slope 2.9 (dashed line). 
In the bottom-left panel both galaxy sub-samples are shown with the best-fit relation forced to have a intermediate slope of 3.20. The 
bottom right figure shows the best-fit Kormendy relation (slope=3.21) produced by brightening the surface-brightness of the z ~ 0.2 
galaxies by 0.6 magnitudes. 



this inherent degeneracy the Kormendy relations can firmly 
exclude very strong growth (> 50% ) from mergers in the 
redshift range 0.2 < z < 0.8. 



6 CONCLUSION 

The results from a thorough re-examination of the BLR HST 
images of a sample of z ~ 1 SCR radio galaxies have been 
presented. It has been shown that, in terms of scalelength, 
absolute magnitudes and Kormendy relation there are no 
significant differences between z ~ 0.8 and z ~ 0.2 SCR ra- 
dio galaxies. The two populations appear to be fully consis- 
tent with being comprised of old stellar populations formed 
at high redshift and evolving passively thereafter. 

It is obviously true that the fact that both the high and 
low-redshift radio galaxies have absolute luminosities con- 
sistent with pure passive evolution can easily be reconciled 
with the involvement of some dynamic evolution. In the dy- 
namical model this simply requires that the proto-galactic 
clumps which merge to produce the final galaxies where 
formed reasonably coevally. The results presented here do 
not then require that these radio galaxies must have formed 



in a monolithic collapse at a single redshift. However, they 
do suggest strongly that the vast majority of merger activity 
within this population of massive ellipticals must have been 
completed before z ~ 1. Although semi-analytical galaxy 
formation models in a f2o = O.S, A = 0.7 cosmology predict 
as much as 70% of present day ellipticals to already be in 
place by 2 ~ 1 (Kauffman & Chariot 1999), it is worth re- 
membering that the same is not true of BCGs, which are 
still predicted to grow by a factor of 2 —> 4 between z = 1 
and the present day (Aragon-Salamanca et al. 1998). 

The crucial question still to be answered is whether or 
not it can be proven that a significant difference in cluster 
environment does exist between the two populations. If it 
can be shown that this is definitely the case then the argu- 
ment forwarded by BLR that we are observing the effects 
of dynamical evolution producing the characteristic mass 
required for powerful radio emission at the two different 
epochs remains tenable; the correlation between luminos- 
ity and scalelength would actually predict the z ~ 0.8 and 
z ~ 0.2 sub-samples to have similar scalelength distribu- 
tions. However, as was discussed in Section hi, while there 
have been numerous surveys carried out recently tackling 
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the environments of high-z radio sources, the work on the 
low-z environments looks to be subject to possible system- 
atic error. This suspicion is further strengthened by our new 
deep R-band HST images of the 10 2 ~ 0.2 radio galaxy sub- 
sample (McLure et al. 1999a, Dunlop et al. 2000). These 
images show numbers of apparent companion objects which 
appear consistent with Abell classes 0^2. This is in good 
agreement with our recent work to model the first to fourth- 
ranked brightest cluster galaxies (BCG) in clusters spanning 
Abell classes 1^4, imaged with HST (Dunlop et al. 1999). 
This work has shown that the scalelengths of these objects 
are not significantly different from the radio galaxy scale- 
length results presented here. 

In an attempt to disentangle the radio galaxy environ- 
ment problem, more near-infrared observations of the low-z 
radio galaxies and quasars from our HST programme are 
planned, to compliment our existing R and B-band images 
of this sample. This combined dataset will allow a fully qual- 
itative investigation of the cluster environments of the low-2 
radio galaxies with which to compare with the existing high- 
z environmental data. 

If it transpires that there is no significant difference in 
SCR environments at high and low-z, then a much sim- 
pler picture of powerful radio galaxy evolution emerges. 
Combined with the scalelength evidence presented here, 
this would suggest that the host galaxies of powerful radio 
sources are basically the same sort of objects at all redshifts 
from 2 = 0^ 1. In this picture the fall-off in SCR radio 
power over this redshift range would simply be due to the 
progressively smaller amount of gas available to feed the 
central engine. 
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